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ABSTRACT 

With a current-splitting technique, we calculate the gauge-invariant initial-state radiation to 
order 0(a) with soft-photon exponentiation for on- and off-shell IT-pair production. This 
result generalizes the convolution formula, which is known from the description of the Z 
resonance, to the case of the production of two IP-bosons. After up to eightfold analytical 
integrations, a sufficiently smooth integral over three invariant masses remains to be treated 
numerically. Including the Coulomb singularity, the largest corrections are covered. We discuss 
the corrections in a large energy range up to -\/^=l TeV and draw numerical conclusions on 
their influence on the hP-mass determination at LEP 200. 

Table 1 and hgure 5 are revised after correcting the analytical formulae for the nonuniversal 
corrections. 
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1 Introduction 


The Standard Model of electroweak interactions allows the calculation of the width Tvi/ [0 and 
mass Mw of the W boson; the latter may be iterated from Ar 0, the electroweak correction to 
the muon decay constant. For hxed t-quark and Higgs masses, the predictions of the Standard 
Model are much more precise than the available experimental data. 

A direct measurement of F^v and Mw will be one of the main tasks of LEP 200 P, |]. This 
will be possible from a study of the reaction 

e+e- ^ (VF+VF-) ^ 4/, (1) 

which is accompanied by the emission of photons and gluons, 

^ 4/ + Tii'f + n2g. (2) 

The photons may be emitted by the initial and intermediate states; and both photons and 
gluons by final state particles. With a centre-of-mass energy at LEP 200 of a/s ~ 190 GeV, the 
VF-pair production proceeds near the threshold. Thus, the finite width effects in the off-shell 
production must be taken into account. And last but not least, virtual electroweak corrections 
have to be inserted properly. 

Here, we will concentrate on the treatment of Initial State Radiation (ISR) 0. The aim is to 
obtain (semi-)analytical expressions for the total cross section at LEP 200, but also at higher 
energies |0, PJ. In section 2, we introduce the Born cross section and define the notations. In 
section 3, photonic bremsstrahlung will be faced. Since the hF-boson is electrically charged, we 
have to derive properly a gauge-invariant definition of ISR. Section 4 contains the cross section 
formulae with contributions from ISR. For the applications at LEP 200, we discuss in section 5 
a possible estimate of the Coulomb singularity. Numerical results and some conclusions both 
for LEP 200 and for a wide energy range are presented in section 6. 


2 The Born cross section 


The main contributions to reaction (|I|) are shown in figure 1: two crayfish diagrams and one 
crab diagram. We use the unitary gauge; in a non-unitary gauge, one should have to take 
into account additional diagrams. The cross section is well described by a twofold convolution 
of a hard-scattering off-shell cross section 0: 





dsip(si) 



ds2p(s2) cro(s; Si,S2), 


( 3 ) 


p{Si) 


1 

TT \Si 


\/G7PvF('Si) 

+ i^iTw{si)\‘^ 


X BR( 2 ), 


( 4 ) 


^In another article we classify and study several background reactions with the same signature as (|^), 
but different intermediate states. They have to be added in order to ensure gauge invariance for the process. 
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Figure 1: The Feynman diagrams for off-shell W-pair production: crayfish and crab. 


The off-shell width r^(sj) contains a sum over all open fermion-decay channels / at energy 
and BR(i) is the corresponding branching fraction. We have rederived the results of [§ 
and refer for details to that article. For later use, we split the cross section into three pieces: 

= o-Q ++ ^0- (6) 


3 Photonic corrections to W-pair production 


In the description of the complicated scattering process (@), it seems to be desirable to hnd 
gauge-invariant subsets of Feynman diagrams. Further, one would like to separate in a reason¬ 
able way the terms related to bremsstrahlung from the genuine electroweak virtual corrections. 

A procedure to disentangle hnal-state real photon emission from the rest of the process (|^) 
has been proposed in l]!^. A corresponding separation of virtual photonic corrections should 
be performed in parallel. In the presence of W exchange, these do not form a gauge-invariant 
subset of diagrams. Being mainly interested in real bremsstrahlung, one could decide to 


combine it only with the singular parts of the virtual diagrams O- The singular parts are 
then uniquely dehned up to an additive constant and will compensate the singularities of the 
real bremsstrahlung. 

These few remarks should indicate that there is much work to be done if one intents to 
perform a systematic, semi-analytical study of reaction (|^). For a Monte-Carlo approach, we 
refer to 0. 

In the following, we will concentrate on ISR. Thus, we avoid to be faced with the full 
complexity of the problems, and at the same time we will cover the bulk of the numerically 
important corrections. The initial-state radiation offers no problems in the case of radiation 
from the crayfish diagrams. Here, the initial charges are separated by the exchange of neutral 
gauge bosons from the rest of the diagrams, and ISR is well-dehned. 


3.1 The t-channel ISR problem — splitting the neutrino current 

The crab diagram contains a neutrino exchange in the t channel. The photonic corrections to 
the initial state are shown in hgure 2. 
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Figure 2: The usual initial-state photonic corrections to the crab diagram. 


The bremsstrahlung may be described by the following electromagnetic current: 


J 'br 

//,a/3 


QeU{k2) 


W2 - k2 

{w2 - k2y'^^ 


—2kip 


2^2 - 
-2k2P 


la 


-wi + ki 

{Wi - 


(1 + 75)m(A;i), 

(7) 


where the {ki, k 2 ,p,wi,W 2 } are the four-momenta of {e“, e+, 7 , hF“, hF’*'}, respectively. This 
current is not conserved, p^J^'^ajS 7 ^ 0) gauge-invariance is violated. Nevertheless, one may 
construct a conserved current also for the crab diagram. Intuitively, it is evident that one has 
to allow the flow of charge within the initial state; a charge Qe has to be associated to the 
neutrino propagator. This, in turn, enables the neutrino to radiate, see hgure 3. The resulting 
auxiliary current. 


/y-aux 


QeU{k2) 


W2 - k2 -Wi ki 

(w2-k2y^^{w,-k,y^^^ 


{l + -f5)u{ki), 


( 8 ) 


ensures current conservation in the initial state and thus makes the d combined diagrams 
gauge-invariant. In fact, the net initial-state current. 


rr'mi — rrhv 1 /T-aux 


(9) 


fulhlls current conservation. 


tJ';! = 0 . ( 10 ) 

In order to compensate for the extra piece in the matrix element, one has to add also a 
diagram with flow of charge —Qe through the neutrino propagator, which in its turn has to 
be combined with the charge flow through the W bosons to build also a continuous electric 
charge flow - but now as part of the intermediate state of the process (or, in case of on-shell 
IF-pair production, of the hnal state). In effect, the electrically neutral neutrino has been 
split into two oppositely flowing, equal charges 0 . 

^For quarks, the auxiliary diagrams from the ‘charged’ neutrino are naturally present, and by no means 
auxiliary. For an application of the CST in ep scattering at HERA, see []ll|. 
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Figure 3: The auxiliary initial-state photonic corrections to the crab diagram. 


For this reason, we call this method of restoring gauge invariance for ISR the Current 
Splitting Technigue (CST). 

The introduction of an additional bremsstrahlung diagram to ISR has a further conse¬ 
quence: one must take into account also additional virtual corrections. These are also shown 
in hgure 3. They, together with the corresponding counter terms, will ensure the compensa¬ 
tion of the infrared divergencies of the photon radiation from the beam particles. The counter 
terms also cancel the UV-divergencies of the newly introduced virtual corrections. The addi¬ 
tional radiation from the neutrino line is free of any infrared divergency, since the radiating 
particle is off-shell. 


4 Initial-state radiation 

For the s-channel crayfish diagram we rederived the ISR correction factor, which is well- 


known from Z-resonance calculations; see e.g. |]^ and references cited therein. The corre¬ 
sponding cross section is described by a threefold convolution: 




ds’ r 


^(H-S')-FFrl p(si)p(s2)cro(s';si,S2). 


( 11 ) 


dsids 2 s 

Here, it is n = 1 — s'/s and the soft plus virtual photon part S{s) and hard part H{s, s') are: 
.2 1 - 


S{s) = 


a 


TT 


Tf 

y ~ 2 


+ H{s, s') — + —j Pe + (12) 


and fde = 2a/7r[ln(s/mg) — 1]. 

4.1 The ISR connected with the crab diagram 

For a treatment of single-photon bremsstrahlung, one has to integrate a tenfold distribution. 
The degrees of freedom of the four-fermion hnal state have been chosen to be two pairs of 
angles in the rest systems of the two fermion pairs. These integrations are trivially performed 
with the method of tensor integration (see e.g. []^). As further integration variables, we took 
the cosine of the photon angle 9.y in the centre-of-mass system and the two angles 0vv, 9w of 
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the virtually produced iy“-boson in the recoil system, i.e. with respect to the photon axis. It 
is the integration over these three angles which demands some effort, and which leads to the 
tremendous complications in the case of the crab diagram with a neutrino exchange in the t 
channel. 

Details of the calculation will be published elsewhere. The algebraic manipulations have 
been performed with the aid of the programs for algebraic manipulations SCHOONSCHIP and 
FORM |T^. We hnally arrived at the following cross section: 


dVa 

dsids2 



ds' 


p{si)p{s2) ^S^ + U, 


(13) 


‘5a(s, s'] Si, S 2 ) 


1 + S{s) 


si, S 2 ) + si, S 2 ), 


Tda.{s, s'] Si, S 2 ) 


H{s, s')cro('S'; Si, S2) + cT^(s, s'; Si, S2). 


(14) 

(15) 


Here, a=st, t denote the st-interference and t-channel contributions, respectively. The func¬ 
tions 5a(s, s'; si, S 2 ) and 7da(s, s'; si, S 2 ) contain extra cross-section pieces with deviations from 
the structure of the s-channel case. As already mentioned, they are nonuniversal in the sense 
that they differ for a=st,t. Further, they are not factorizing, i.e. they do not contain the off- 
shell Born cross section as an explicit factor. A third feature of them is a screening property: 
compared to the Born cross section, they have a damping overall factor. 


st,t / ! \ ' 5 iS 2 

a^^^(s ;si,S2) 


( 16 ) 


which suppresses potential mass singularities, which otherwise would be generated by new 
kinematic logarithms compared to the Born cross section. The non-universality of the ad¬ 
ditional corrections could, in principle, spoil the gauge cancellations 10 between the three 
diagrams of figure 1. From the analytical structnre of the nonuniversal terms, such a can¬ 
cellation may not be seen; the (t| ^{s']Si,S 2 ) are complicated functions, as one may expect 


after seven to eight integrations. The property ( [TT|) is welcome in this respect, but it could be 
compensated at least partly by the integrations over invariant masses; numerically, we have 
observed that the net cross section behaves properly at extremely high energies. A conse¬ 
quence of the auxiliary current is that our formulae for ISR contain some contributions which 
in other calculations are considered a part of the intermediate or final-state radiation. 

Besides rational fnnctions and logarithms of s,s',Si,S 2 , they contain dilogarithms and, in 
case of the virtual corrections, trilogarithms. We will publish the explicit expressions for the 
cross section elsewhere. Despite the fact that they are relatively lengthy, we would like to 
stress that they are much more compact than they would be without the introduction of the 
auxiliary current (^). It is of importance for a numerical handling that the cr| and are 
small in the LEP 200 energy region. This allows us to include soft-photon exponentiation and 
0{a‘^) leading-logarithmic corrections into the st- and t-channel universal terms, and to treat 
the nonnniversal rest to order 0(a). 

At the beginning of this section we reproduced, in the context of off-shell FF-pair produc¬ 
tion, a general resnlt for s-channel ISR: the radiator fnnctions S and H, which describe the 
creation of an intermediate vector boson V* together with a photon, e’''e“ —R*(s') 7 ini. The 
corresponding t-channel contribntion with two virtnal vector bosons. 


e+e- ^ R;(si)R;(s2)7ini, 


(17) 
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and its interference with the s channel, is characterized mainly by the same, nniversal fnnctions 
S and ff, but also by some deviations S and II, due to the more involved kinematics. It was 
explicitly assumed that and V 2 have no common hnal-state interactions; such interactions 
would destroy the anticipated kinematical situation. 


5 The Coulomb singularity 


In the threshold region, there is another type of photonic corrections besides the ISR, which is 
potentially large: the Coulomb singularity []^. For off-shell hF-pair production, it originates 
in the s channel from the insertion of a virtual photon to the and ZW^W~ vertices; 

in the t channel, it is due to the vW^W~'^ box diagram. It is not completely clear to us how 
to take the Coulomb singularity into account numerically, apart from the fact that it should 
appear as a universal, positive factor in the net cross section, which takes into account the 
off-shellness of the process and should approach the known on-shell correction. 


ac — <j[l + C{s'] Si, S 2 , M^r^IIw)\ I ( 18 ) 

C(s'; si, Mw, Fm. ^ 0) ^1^, P{s') = ^Jl-AM^/s'. (19) 

With a Feynman-diagram calculation one may hnd the scalar one-loop function, in which the 
Coulomb singularity is located. Unfortunately, an evaluation of it without a knowledge of 
the kinematic behaviour of the imaginary part of the W propagator is impossible. A crude 
estimate, assuming for the IF-propagators (§) under the loop integral that i^/sirw{si) ~ 
iSiTw/Mw, yields for the leading term in the small-width approximation: 


C{s-, si, S 2 , Mw, Fiy) — [ 1 — 3 


^ w 

Ml-. 


TTa 


2/5(s';si,S2)’ 


/3(s'; si, S 2 ) = - V A(s'; Si, S 2 ). (20) 


Here, A is the usual kinematic A-function. 

A different approach which is based on a non-relativistic treatment of the threshold region 
has been developed in |^, 


For details we refer to the original literature. 


6 Numerical results 


In the discussion of numerical results, we will concentrate on two observables, which may be 
used for a precise determination of the W mass at LEP 200 @]. The calculations have been 
performed with the aid of the FORTRAN program GENTLE []^. The corrected numerical 
values have been produced with version CCll of July 1995. We have used the following input 
parameter dehnitions, which slightly deviate from what was used in 1993: 


G, = 

a = (y(^2Mw) = 


’w 


= s 


1.16639 X 10 
1/128.07, 


-5 


2.e// 

Tia 

W 


Mz 

= 91.1888 GeU, 

Fz 

= 2.49747, 


GeV 


-2 
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( 21 ) 


M’\y 

Tw 


80.230 Gel/, 
9G^M^ 
67ry2 


We will discuss the radiative corrections both for the cases of on-shell and off-shell W- 
pair production. Our ‘minimal’ radiative corrections are the universal corrections introduced 
in (p!2[), but neglecting the indicated terms of the order of C>(q(^) [Q. Additionally, we may 
take into account: 


(i) the universal 0{a^) terms in ([T^) [f/ 2 = 0 ,l]; 

(ii) the nonuniversal cross sections and afj in (|T^ and (pAl), which are the main result of 
this study [^=0,1; .^=0,1]; 

(iii) the estimate of the Coulomb singularity [G=0,1]. 


6.1 The IT-pair excitation curve 

We should like to start the discussion with a few comments on figure where we show the 
cross sections crg“ and for Born on- and off-shell W-pair production, and also for 
the off-shell production with universal ISR (G 2 =l, S=H=C=0). Compared to the cross 
section develops a tail. Although this tail becomes heavily suppressed at high energies, it 
is not vanishing i. The other tail phenomenon is due to the universal ISR and is much more 
pronounced, but again weaker than e.g. is observed for the narrow Z-resonance shape. In the 
hgure, the radiatively corrected on-shell cross section is not shown. As a matter of fact, 
we mention that the relative differences between Ug'^ and cxg®^, and between and are 
quite similar. 

Figure ^ and table allow to estimate the differences between the simple universal ISR 
approximation and the exact treatment of ISR, which additionally includes the nonuniversal 
cross-section contributions connected with the crab diagram, and also the Coulomb singular¬ 
ity. The nonuniversal corrections and afj- increase the total cross section by -|-0.4%; -|-0.8%; 
-|-1.5% at a/s = 165-i-190; 500; 1000 GeV, respectively. We quote these values for the on-shell 
case; the off-shell corrections (in percent) are nearly the same. 

On top of that, the Coulomb singularity yields a positive correction, which has its maximum 
value of about 6% at the threshold and asymptotically approaches the on-mass-shell value of 
^Tia = 1.15%. At large a/s, in contrast to the threshold region, it is only one of many 0(a) 
corrections. We compared numerically our equation (^) with the predictions of equations (1) 
and (9) i of [|T^], and got very good agreement in the threshold region. The predictions of 
the non-relativistic calculation are slightly smaller, but the absolute differences of the two 
calculations do not exceed 0.5%. 


6.2 The radiative energy loss 

Two potential methods for a determination of Myy are the direct reconstruction of events and 
the measurement of the (upper and lower) energy end points in leptonic IF-boson decays. 


^In [^, with which we analytically agree, this is not seen from the figures. 
^In equation (9) of |Q, we treated the arctg as defined between 0 and tt. 







Figure 4: Total cross section for W-pair production: as and universal ISR corrections. 



Figure 5: Total cross section for W-pair production with ISR and Coulomb corrections. 
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H 

U 2 

c 

a, pb 

(Erad), GeV 

0 

0 

0 

0 

16.405 

2.091 

13.634 

1.158 

0 

0 

1 

0 

16.401 

2.153 

13.602 

1.190 

1 

0 

0 

0 

16.471 

2.091 

13.686 

1.159 

0 

1 

0 

0 

16.405 

2.090 

13.634 

1.158 

1 

1 

0 

0 

16.470 

2.091 

13.686 

1.159 

0 

0 

0 

1 

16.783 

2.107 

14.042 

1.170 

1 

1 

1 

1 

16.846 

2.170 

14.063 

1.201 


Table 1: The off-shell W-pair produetion cross section and the radiative energy loss at two 
energies: i?beam='95 GeV (upper rows), 88 GeV (lower rows). The S and H switch on and 
off the soft and hard parts of the non-universal QED corrections, the U 2 the second order 
corrections, and C the Goulomb correction. 
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They both rely on the knowledge of the effective beam energy, which deviates from the beam 
energy itself by the radiative energy loss -Erad- Its average value is 


{Erad) 


ko 




dsids2 


dfa 


ds'dsids2 ’ 


( 22 ) 

(23) 


Since the radiative energy loss is essentially due to ISR, with a possible smaller amount from 
initial-hnal interferences, one may assume that the formulae of the foregoing sections cover 
the bulk of effects. 

Numerical estimates for the cross section and the radiative energy loss are collected in 
table 1 for two typical energies at LEP 200. The radiative energy loss amounts to about 
2200 MeV at y^= 190 GeV (1200 MeV at ^^=176 GeV). These shifts are mainly due to 
the universal, 0{a) ISR corrections with soft-photon exponentiation. Further, there are at 
190 GeV shifts of -|-62; ±0; -|-16 MeV due to: universal, non-exponentiated higher orders; 
nonuniversal corrections; the Goulomb singularity (at 176 GeV: -|-32; ±0; -|-12 MeV, corre¬ 
spondingly). Apparently, these corrections have to be taken into account properly, if one aims 
at an accuracy of 50 to 100 MeV for the fV-boson mass. 


The numerical results of were based on the corrected version of the Fortran program. 


6.3 Conclusions 

We have performed the hrst complete, gauge-invariant calculation of initial-state radiation cor¬ 
rections to off-shell IV-pair production. In the threshold region, we also include the Goulomb 
singularity. The net improvements compared to the well-known universal ISR corrections are 
relatively small, but non-negligible, at LEP 200. They become more and more important 
at higher energies, and the hnal-state radiation seems to be much less pronounced after the 
inclusion of the auxiliary current into the ISR. 
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